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1 Introduction 

Consider the Wigner random matrix ensemble X^v = where 



d \ 2. Xii are iid real-valued with mean and variance 2/N, and 



1. Xij, i < j, are iid complex v alued with mean and variance 1/N. 



for i < j, i.e., the matrix is Hermitian. 

, 4. The distributions of Xij and Xu are supported on at least three points, and have finite 

O ! fourth moments. 



The main result of the present paper is 

Theorem 1.1. Let fi, . . . , fp be compactly supported functions on R*^i, . . . , R'^p respectively, 
—2 < ai < . . . < ap < 2, and li = [a^ — e^, + e^] such that li 's are disjoint from each other. 
For a compactly function g on M.^ , and a G (—2, 2), denote 



: g[a]{Xi, . . . , Xn) = Yl g{Npsc{a){Xi,-a),...,Np^c{a){Xi^-a)). 

l<n<...<ifc<Af 

Then 

limEjJ— / /i[Mi](Ai,...,AAr)ciMi- JJE— / fi[ui]{Xi, . . . ,XN)dui = ^. 



Let Ai := A\^^ = [Ei, Ei + j^], i = 1, p for any fixed p, where Ei are distinct points 
in (—2, 2) and q are fixed real numbers independent of N. Also let A/'(A) be the number of 
eigenvalues in A. 

As a corollary we have 
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Corollary 1.2. The following averaged asymptotic independence of eigenvalue counting 
functions hold (see the Theorem 3 of IBiDeNaf for GUE) 

lim ^ [ ... f P[A/'(Ai) = /i, . . . ,U{\) = l,]dE, ...dE, 

-Utt m{A,) = k]dE, = 0. 

Vast amount of efforts has been devoted to proving that the eigenvalue distributions of 
Xtv as — )■ oo have the same limit as that of the Gaussian Unitary Ensemble (GUE), where 
d of a particular distribution, namely Gaussians. Whenever a property P of the 
GUE limiting spectrum is also true for any Wigner ensemble defined above, we say that P 
is universal. 

For single eigenvalue marginals, one can prove that the empirical spectral measure YliLi 
converges almost surely to the semi-circle distribution, whose density is given by psc{x) = 
^-\/4 — x'^l\x\<2- Thus the eigenvalues of Xjq are supported on the interval [—2,2] in the 
limit. 

The next order of investigation is the joint distribution of more than one eigenvalues. In 
the case of Hermitian ensembles, one could label the eigenvalues from 1 to according to 
their relative orders on M. This perspective is very useful in proving the universality of local 
eigenvalue distribution, to be defined more precisely below (see |ERSYj ). 
A more natural point of view is to count how many eigenvalues lie in a region A C [—2, 2], 
the size of which could shrink as A^ — )■ oo, in order to have a nondegenerate limiting distribu- 
tion. More generally, we define the so-called m-point correlation function Pm{xi, . . . ,Xm) = 
pl^\xi, ...,Xm)hy 

/ f{Xi,...,Xni)Pm{Xl,...,Xm)dXi...dXm = 'K ^ / ( Aii , . . . , Aj^ ) , 

l<ii<...<im<Af 

where / is any bounded continuous test function on MJ^. Notice that pm is not a probability 
density on since it integrates to (^j^^'^^y against the constant 1 function. 
It turns out that for GUE as well as a large class of unitarily invariant Hermitian en- 
sembles(but not Wigner matrices in general), these m-point correlation functions can be 
expressed as a determinant of an m-dimensional matrix, whose elements are given by a fixed 
kernel as follows: 

p^UE(N)(^^^ . . . , x^) = detiK^ixi, x,))™=i, 
where Kj,ix,y) = N ^^i-^)^N-^iYVN)-^^.,i.VN)^r.{y^) ^ ^ ^l£lLhu{x) is the Her- 

mite polynomial /i^ normalized and multiplied by square root of the Gaussian weight. 
Thus in order to study the limiting behavior of as A^ — > ooj it suffices to understand 
KN{x,y). It turns out, using asymptotics from orthogonal polynomial theory, 

1 , c , sine 
lim — -——K]sf[x, X + —j = . 
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In the past decade or so, through the various effort of Brezin and Hikami |BrHij . Johansson 
[JoOlj . Erdos, Ramirez, Peche, Schlein, Yau, Yin |ERS Yj . |EYYj . and Tao,Vu |Tao Vuj . 
the universahty of the sine kernel hmit in the bulk of the eigenvalue spectrum have been 
established for very general Wigner ensembles, with only a mild condition of subexponential 
decay for individual entry distribution. 

In this article, we look at a different property of the eigenvalues, namely independence in 
the limit — )■ oo of two or more eigenvalue clusters located at macroscopically separate 
positions on (—2,2). According to our best knowledge, the only results in this drection so 
far are given by 

1. Bianchi, Debbah, and Najim |BiDeNa] . who showed that 

p 

lim |P^[Ar(Ai) = h,... ,Af{Ap) = Ip - n P^,[Ar(A,) = /,]!= 0, 

N—>-oo -*- 

1=1 

where Aj = [Ei,Ei + ^] for £"1 < E2 . . . < Ep G (—2,2), i.e., distinct locations 
in the limit spectral support, and q G M. If q are negative, then we interpret 
Ai = [Ei — ^, Ei\. Finally A/'(A) denotes the number of eigenvalues in A, here under 
the GUE ensemble only. In fact they showed more: if Ei = —2 and/or Ep = 2, then 
the same result above holds provided Ai = [—2, —2 + jfr/s] ^p = [2, 2 + -^t/s]^ i^i 
accordance with the scaling at the spectral edge. 

In the appendix, we reprove their result in the bulk using the inclusion-exclusion prin- 
ciple, instead of their more sophisticated approach using Fredholm determinant and 
complex analytic tools. 

2. F. Boremann |BolOj reproved the asymptotic independence of the extreme eigenvalues 
Ai and \n for GUE using abstract operator theoretic argument. He also obtained 
the next lowest order term in their joint distribution, and used it as a distinguishing 
statistics for GUE against other Wigner ensembles. 

3. Soshnikov's earlier result |So99] on universality of the edge distribution of Wigner 
matrices whose entries have symmetric distribution and satisfy a sub-Gaussian moment 
assumption actually proves the asymptotic independence of Ai and Aat as well, using 
his high power moment method. 

4. Borodin, Okounkov, and Olshanski [BoOkOl] established the asymptotic connection 
between the bulk of GUE spectrum and the Frobenius coordinates of a random parti- 
tion distributed according to the Plancherel measure. They mentioned that the distri- 
bution of Frobenius coordinates at far away regions are asymptotically independent, 
which presumably is proved in their paper also. 

In the present article we show that the result from [BiDeNa] is universal. 

2 Outline of the proof 

We borrow heavily ideas from |ERSY] . We illustrate the proof with p = 2 only. The more 
general case follows the same pattern, without additional assumptions, unlike the main 
results in [EPRSYj . 
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The first step is to show that the eigenvalues confined in two clusters Ai U A2 := 
{^Li+i, ■ ■ ■ , ^Li+n} U {\l2+i, ■ ■ ■ , Aij+n} couvcrges in total variation distance to their joint 
stationary measure conditioned on the positions of the remaining Aj's A \ (Ai U A2), under 
the mean- reverting Dyson's Brownian motion in time t << 1. Here n = for some small 
5 E (0, 1), and < limTv < limN < 1, so the two clusters are anchored at distinct points 
on (—2, 2) in the limit. 

As in |ERSYj . we prove convergence in total variations using relative entropy bound by 
establishing a log-sobolev inequality for the invariant measure, which can be achieved very 
elegantly using Bakry-Emery theory. We will prove that the Hessian of the corresponding 
local Hamiltonian has spectrum bounded below by a constant in the next section. 
Next we show that for almost all Li < L2, the conditional equilibrium measure of Ai U A2 
is close in total variation to a product measure on Ai x A2, where the components are the 
invariant measure on Ai conditioned on A \ (Ai U A2) and that on A2 conditioned on the 
same event, for asymptotically almost all A \ (Ai U A2) configurations. 
In [ERSY] . it has been proved that the correlation kernel of Ai (or A2) when rescaled so 
that suppAi = [—1, 1] converges locally to the sine-kernel. This relied on results from |LeLu] 
where it was shown that sine-kernel is universal in a broad class of unitarily invariant en- 
sembles whose matrix element joint density is given by exp(— A^TrH(X)). 

Notice that Li and L2 are the ordinal index of the eigenvalues, which are not directly 
related to their rescaled positions in [—2,2]. However, a concentration result gives that the 
eigenvalues do not deviate from their semi-circle law predicted positions by more than 0{N). 
More recently, [EYY] has shown that the deviation cannot exceed a polylogarithmic amount 
with high probability uniformly over all eigenvalues up to the edge. So one can essentially 
treat Li, L2 the same as Ei, E2 the corresponding energy levels. 

Finally under the framework of [ERSY] as well as their sequels, an averaging over a 
macroscopic range of energy levels Ei and E2 is needed for the above sine-kernel convergence, 
because the convergence to local equilibrium cannot be shown for all Li. 

Thus their result was of the form 



lim lim — / dE dadbO(a,b) , , , ^ P2(-^H — ,E , 



(1) 



g{u)[l - {^^r]du, (2) 

TTU 

where 0(a, b) := g{a — b)h{^^) and g, h are both compactly supported with J h{x)dx = 1. 
Note that h is just a local bump function that does not appear on the right hand side. If we 
do have 

1 , ^ a ^ b , .smnu.n 

;P2{E+ E + — -) = 1 ^ ^2 



then ([T]) would follow by a change of variable (a, b) i-> (a — 6, a + 6). 

Using the very recent rigidity result of Erdos, Yau and Yin |EYY] , one gets that the support 
of Ai, [Xli, Ali+„+i] contains a subinterval of [—2, 2] of size N''^ for some 7 < 1 — ^ (recall 
n = N^) with high probability. Hence provided we know Li is a good index, sine-kernel con- 
vergence without averaging can be established in the interval [7^-^(Ll) + ^N^-^ ,n-^{Li + 



4 



n + 1) - IN^-^]. 

Moreover, if we have two eigenvalue clusters Ai and A2 based at two good eigenvalue in- 
dices Li,L2, they would contain intervals Ii := [Ei,Ei + N''^] and I2 := [£'2,-^2 + N~^] 
respectively with high probability. So conditioned on A \ (Ai U A2), which stays away from 
Ji U I2 with high probability for the same reason, the distribution of the counting functions 
l_V(Ai)=h and ljV{A2)=h with Ai C /i and A2 C I2 are almost independent for almost all 
A \ (Ai U A2) configurations, by local convergence to equilibrium. 

This is however not yet enough to show lj\f(Ai)=h W{A2)=Z2 independent uncondi- 
tionally. For that one needs 

F[^^{A^) = /i, Ar(A2) = /2IA \ (Ai U A2)] = fih, h) + 0(1) 

asymptotically almost surely, for a fixed function /. The relation between universality 
o£A/'(A) and that of the local correlations, as established in the appendix, reduces the prob- 
lem to showing that pm\xi, . . . , a:rn|A\ (Ai UA2)) for Xi G Ai and for Xi G A2 both converge 
to a constant function almost surely for each m. But this is essentially proved in |ERSY] 
and strengthened by results in |EYY] mentioned above. Finally we integrate over all config- 
urations of A \ (Ai U A2), using the boundedness of the integrand, to obtain 

E[P[Ar(Ai) = k,Af{A2) = /2IA \ (Ai U A2)]] = fih, I2) + 0(1) 

This establishes the universality of asymptotic independence for ERSY-ensemble (i.e., Jo- 
hansson matrix with vanishing GUE component). The crucial assumption for the above 
argument to work is the 3 conditions (2.4)-(2.6) imposed on the entry distribution of the 
Wigner component of the ERSY-ensemble. Those ensure that the global entropy is bounded 
by A^^^*^ for e > 1/4 at the end of a burn- in period of t = 0(1/A^). 

The next step is to remove the Gaussian divisibility condition by invoking the powerful four- 
moment theorem of Tao and Vu, which states that if two Wigner ensembles and have 
entries with the same first four moments, then their joint eigenvalue distributions for finitely 
many eigenvalues are very close in the sense of weak convergence. 

Theorem 2.1. ( ITaoVi^ ) Given that and Yat match up to fourth moments, and let Xi 
be the eigenvalues of Xn and those ofY^, then for some universal constant K,c> 0, 

|EG(A.„ . . . , A.J - EG(A:^, . . . , A:J| < Kn-^ 

for all G with V^G{x)\ < n^^'^, where j = 1, . . . , 5. 

Note that the version above differs from the original theorem of Tao and Vu by the power 
of n in the derivative bound of G because they scale Wigner matrices by a factor of N. 
If we take /(xi, . . . , X/J supported on Ai and g{xi, . . . ,xi^) supported on A2 > Ai, satisfying 
the derivative condition above, and denote 

nfg) ■■= ^'^' e Y1 /(^-' • • • ' KM^n, ■ ■ ■ , A, J, 

ii<...<il-^<ji<...<jl^ 

then 

\ExM9)-^yAf9)\<KN-'^. 
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From here it is easy to derive that the correlation function of the form 



(iV-/i-/2)! 



iN-h)\ 



m 



for Xi G Ai and yj G A2, provided same holds for Xj^. Thus if a matrix ensemble has the 
same first four moments as some Johansson matrix ensemble, it also enjoy the property of 
asymptotic independence. As pointed out in |Tao Vu] . if the entry distribution is supported 
on at least 3 points, then one can always find a Johansson matrix that match the first four 
moments, which concludes the proof. 

3 Convergence to local equilibrium for more than one 
eigenvalue clusters 

In this section, we prove convergence in total variation distance of the joint distribution of 
two mesoscopic eigenvalue clusters to the conditional equilbrium measure given the other 
eigenvalues. For ease of notation, we will stick with two eigenvalue clusters in this section. 
Our notation will be largely based on |ERSY] . except that we will denote by := (1 — 
ty/'^W + t^/^f/ the Johansson matrix ensemble. Let G [N] be two sequences of 

integers such that limTv = a and limjy -^q- = b both exist and 



Let n = where 6 G (0, 1) is fixed, but can be chosen as small as one likes. For simplicity, 
we will suppress the dependence on N in the above notation whenever possible. Consider 
the following two contiguous eigenvalue clusters 



Clearly these two clusters are disjoint for all large N. Define the relative entropy between 
two probability measures on the same space by 



whenever / exists, where / = ^ is the Radon-Nikodym derivative of u with respect to /x, 
and S{fi, u) = 00 otherwise. 

Also define total variation distance between /i and u by 



< a < 6 < L 



A2 := (Al2 + 1, • • • , Aij+n)- 




fi — z/||tv = sup \fi{A) — i^{A)\. 

Acn 



When / exists, one can alternatively define 
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Thus one expects total variation distance to be bounded by relative entropy. Indeed we have 



For the reader's convenience, we include a proof of this fact in the appendix. 

It turns out that entropy analysis is well-suited for Dyson Brownian motion, via the so- 
called Bakry-Emery theory, which leads to a logarithmic sobolev inequality that relates the 
relative entropy with Dirichlet form applied to the square root of the running density of the 
Markov process (see |GuZej for an excellent introduction to the subject; also [ERSY] section 
5 for its application to DBM). 

In [ERSYj three conditions (2.4) - (2.6) are assumed about the initial Wigner ensemble, 
which basically stipulates that the distribution of the individual matrix entries satisfies a 
logarithmic sobolev inequality and has sub-exponential decay, and that its fourier transform 
has polynomial decay for a sufficiently large polynomial degree. 

These distributions are flexible enough that their flrst four moments can be anything that 
comes from a distribution with at least three points in the support (see [TaoVuj Corollary 
30). We will assume the same set of conditions in the next three sections. Under these 
conditions, we have the following bound on the global entropy after a burn-in period of 
0{1/N) of the eigenvalues of the Wigner copy W evolving under the DBM: 

Lemma 3.1. (Lemma 5.1 of lERSY^ ) Let ft he the density of the eigenvalue distribution 
of Xj^ := (1 — ty/'^W + t^/'^U with respect to the GUE spectral measure fx{dXi, . . . , dX^) = 
Z^^e"^^'^i^?"^'<^^°^(-**'"^^^'dAi . . . dXN, and where U is a standard GUE{N) matrix. Then 
for every a > 1/4, 



Next using (5.4) of [ERSYj . we can bound the global Dirichlet form of time 

t > N-^ by 



for t = tN-^ and r > 2 (see (6.1) of [ERSYj ). 

The reason to work with Dirichlet forms instead of entropy directly is that the former has 
convenient additive property, which allows us to spread the global Dirichlet form into local 
ones. Let A denote the vector of all eigenvalues of a Hermitian matrix ensemble, arranged 
in increasing order. The notation A \ (Ai U A2) thus has obvious meaning, i.e., the vector of 
the external eigenvalues arranged in increasing order, and so does A\ Ai. We will always use 
the letter y to stand for a realization of A \ (Ai U A2) and x^^-* and x^"^^ to denote that of Ai 
and A2 respectively. Also we denote collectively x = (x^^-*, x*-^-*). Unlike in [ERSY] . we index 
the components of y consecutively. In analogy with the deflnition of localized Dirichlet form 
in section 6.2 of [ERSY] . we deflne the following 2-point conditional local Dirichlet form by 



~ ^IItv — 25'(/i, z/). 



S{fN-ifi,fi)<CN 




This easily gives 




,ii,i2 2A^ 



V./|W^'^^(x), 
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where Hy^'^^ is the invariant measure on Ai U A2 given A \ (Ai U A2) = y and 

where for a subset J C R, we denote by (/")^ the set of increasing n-tuples with components 
in /. Note that [yL2-n, yL2-n+i] is the confining range of A2. Denote by ft,y{x) the conditional 
density of the eigenvalues A1UA2 at time t given the external configuration y, with respect to 
the conditional invariant distribution on Ai U A2 given y. By exactly the same computation 
as in [ERSYj between (6.17) and (6.18), we have 

w n f /VT~\\ ^ V" /"^ ^ \^xjtiy,x)\^ + \^x,+„My,x)\^ 

Ef,D^L,,L2),y{\l ky{x)) = ^ l^J dxdy 'U^) u{y,x). (3) 

where u{y,x) is the GUE eigenvalue density with respect to the Lebesgue measure on the 
positive orthant in M^. The expectation above is taken with respect to the external config- 
uration y. 

Now we sum the above equation ([3]) over Li G [Na — n2, A^a + n2] and L2 G [A^6 — ^2, Nb + n2] 
where n2 = N^^ and 5 < ^2 < 1. This amounts to summing the Dirichlet form for each indi- 
vidual particle in the range at most n times. Thus we get 



^ L\<^[Na-n2,Na+n2\ L2<^&[Nb-n2,Nb+n2\ 

for some constant C. Thus by the discrete Markov inequality, if we let 
Q^{a, b) := {(Li, L2) G [Na - ^2, Na + ns] x [Nb - ^2, Nb + ^2] : Ef^D(^L^^L2),y{\/ kyi^)) 

be the set of good pairs of eigenvalue indices, then for some constant c. 



\gN{a,b)\ > ^ _ c 
4^2 ~ n 



Furthermore, if Fi,,^, := {y G (M^-2")t ; D^L,M),yi^y) < CN^+^^j-'}, then by another 
apphcation of Markov inequality, for (^1,^2) G ^Ar(a, 6), 



Next we compute the conditional local Hamiltonian T-ty{x), which defines the conditional 
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local equilibrium measure iJ>y{dx) — exp{—'Hy{x))dx: 



i=l l<i<j<n l<i<j<n l<i<j<n 



N-2n n 



(4) 



ie[n]je[n] k=l i=l 



We can compute the Hessian of Hy as follows: 



(5) 



n N—2n ^ 



d^.Uyix) = ivxW - E (1) _ . (1) - E ji) I - E (1) _ • 

ie[n]\{i} -^i "^j j=l -^i -^j k=l -^i Uk 

Hence 



+ E /^(l) _ ^(2)^2 ^ ^ TJI) 

2 



d^(2)d^mny{x) = 



The other mixed partials d^(2)d^(i) and d^(2)d^(2) are obtained similarly. To obtain a spectral 
lower bound for HessT^y, we take a test vector v e M^"', 

n 1^2 N—2n ^ 

{V, Ressny{x)v) = Et^ + E . a)_^(i)p + E . (1) _ ^(2)., + E . (i)_^.2 ]^- 

i=l ki^i -^k ) k=\ V'i j I k=l l-^i ilk) 

n -, n ^ N—2n 



2 



+ Et^ + E . (2) _ (2).2 + E . (2) _ (l)^2 + E . (2) _ ^2^^^^+^^ 
i=l jfcT^i V-^i -^ik ^ k=\ -^j ) k=l l-^i yk) 

n n 

^ 2^ _ J1)n2 + / (2) _ (2).2J ~ ^ 2^ _ ™(2)n2 

«=1 j^i V'^'i •^j J {•^i •^j J i,j=l {•^i •^j J 

= m\vf + f^czi J""' + (25" ) + E E[ ar'^m + 52r''iy^' ] 



+2E 



Ji) _ Jl)^2■ 



Thus for X = with x(^) e ([t/Li, yLi+i]")^ and x^^^ e 

N-2n 

Hess?^v(x) > N + inf C V 7^ n2 • 

— myLo — n 
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By rigidity of the eigenvalues of the Wigner ensemble Xj^, proved in |EYYj . the set Q2{Li, L2) : 
{y G (R^^^n^t : max{?/ij+i - y^^, y^^.^+i - i/La-n} > ^} has vanishing probability. Indeed 
Pxjv(^2(-^i, -^2)) = 0(iV^^), for some e > 0. From this we easily conclude 

Hess?^v(x) > — —, 

for almost all y. 

Next we invoke the Bakry-Emery theory as in |ER,SYj on the generator of the DBM 
£ := 2^(A - {VHy, V)), which gives 

dtD^L^,L2),yWft) < —^D{^/ft), 

where again ft = e~*^/o is the the global density of the DBM at time t, and /o is the 
density of the joint eigenvalue distribution of starting ensemble Wn- This further implies 
the following logarithmic Sobolev inequality for any probability density g 

CN 

D{^9) > -^S{g), 

where we abbreviated S{g) := S{g, fiy^'^'^ , fiy^'^^). See Appendix B for a proof of this. Note 
that the logarithmic Sobolev inequality bears no direct relation with any dynamical process 
that converges to the measure fiy^'^^- In fact we don't want to treat it from a dynamical 
point of view since the environment y is changing during the time interval [0, t] under DBM. 
Thus since D^LMyiVTt) < CN^+"^t-\ , we get for {L,, L^) E GnM, 



„6 

S{ft) < CN'+'^^T-' 



So if r >> N^'^'^n^ /n^, which is still o(l) if we choose n2 sufficiently large and n sufficiently 
small, we get S{ft,yfiy^'^^ , fi^^'^^) = o(l), and hence 



4 Asymptotic conditional independence of distant eigen- 
value clusters 

In this section, we show that the invariant measure /^^^'^^ asymptotically converges to the 
product measure /i^ ^ ® fiy^ with high probability. More precisely. 

Proposition 4.1. For each N there is asetfl'^''^^ C ([-2, 2]^-2")t with F[A\{A^^ UA^^) e 
Qjy] —> 1 such that for (Li, L2) G Qn{cl, h) the set of good index pairs, 

for some e > and universal constant C . Here z^^^ G ([^/Li, ^/Li+i])^ and z^"^^ G {[yL2-n,yL2-n+i]) 
are arbitrary. 
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Proof. Recall the 2-point local equilibrium measure anchored at Li,L2 has the following 
y-conditional density 

where l-iy^x) := Tiy^'^^^x) is given by (j4]). The key observation is that 'Hy^'^^{x) splits 
approximately into a sum of two decoupled Hamiltonians: 



where 



n ^ N—2n n 

i=l ^<i<j<n k=l 1=1 



^ n n 



and 



fc=l i=l 
" 1 1 

i=l l<i<j<n 
^ N—2n n ^ n n 

k=l 1=1 k=l i=l 



n n 

2 



^^(^) = E E ^°g(^* ^ ^fc^^)^ + iog(x,+„ - 4^^)^ - E E ^°g(^^ ~ ^•?+") 

A;=l i=l j=l j=l 

First observation is that Rz{x) has small fluctuation for almost all y, i.e., 

sup R^ix) - inf < sup C E + ""^^ (1) ^ (6) 

^•'^ fc=i i=i - 4 3;i+„ -zl' 

where we used the first derivative estimate on the function log(xj — ZkY, and the rigidity 
result from |EYY] on the external eigenvalue vector y to lower bound the separation Xi — z^jf' 
etc. 

Next let Zi, Z2 be the partition function for "^^^^(2) and "H^^^jij respectively, i.e.. 



'([yi^,yi^ + i]")t 

Then comparing the following two quantities 



Zi= exp(-H^^^(2)(a;))c/a;. 



ZIZ2 '■ 

z ■= z 



I exp(-?/^^^(,)(x) - n';^^,,,{x))dx^'^dx^^^ 



N,{Li,L2),y 



[ [ exp(-H^^^(2)(a:) - n^',,,ix) + R.{x))dx^'^dx^'\ 
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we find for a fixed z vector and any x that 

ZiZ2exp{R^{x)) exp{-C—) < Z < Z1Z2 exp(i?^(x)) exp(C— ). 

This imphes that the ratio of the density of /x^ ^'^^ with /^^^^{2) ® l-^yuzW satisfies the following 
bound for almost all z, 

, Cn^^ Z7^exp{-'Hi{x^^^))Z^^exp{-'H2{x^^^)) ,Cn\ 
^"P^-l^^ ^ Z-^expi-nix)) ^ ^"P^l^^- 

Using now the interpretation of total variation distance, we easily arrive at 

/Cn^ 

< exp(^) - 1. 

□ 

5 Unconditional independence 

In this section we will fix a sequence of good pairs of indices (Li, L2) E Q]\f{a,b). By applying 
the main result from |LeLu] . |ERSY] showed that the process on x*^^^ = {x^l^ < . . . < Xn^) 
under the measure fiy^{dx) = exp{—'Hy^{x))dx has the bulk correlation that converges to the 
sine-kernel when appropriately scaled so that = —2 and i/Li+i = 2, for all y in some set 
Y^^ that has asymptotic probability 1 under the spectral distribution of the ERSY ensemble 
= (1 - ty/^WN + t^^'^HM defined before. Here Hy^ix) is defined as before by 

^ n ^ ^ N—n n 

Uy{x) = N[-^x^i - — ^\og{xi- XjY - — J2^^og{xi - VkY]. 

i=l i<j k=l i=l 

Similarly one get sine kernel for x^^^ = {xf'\ . . . , Xn^), for y G Y^'^ with asymptotic proba- 
bility 1. 

The reason is that the Hamiltonian above is essentially the trace of some analytic function 
of the matrix ensemble, hence falls in the category of unitarily invariant ensembles, whose 
local correlation has long been conjectured to converge to the sine kernel, and proved under 
various conditions. To be explicit, we define the rescaling map 



^ yL2-n+l — yL2-n 



Furthermore let Ty^\x''^^) = {tI}^\x^^^), . . . ,Ty^\xn'^)) and similarly define Ty^^(x^^''). Then 
we have 
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Lemma 5.1. There exist sets Y^^'^^ C ^" with asymptotic probability 1 such that G 
yj^^'^^'^ implies sine-kernel universality of local correlation for both Ty^\x^^'') and Ty'^\x^'^'^). 
More precisely, for a compactly supported continuous function f , and a G (—2,2), 

^VN XI fi.'^pM{Ty{xf^) - a), . . . ,np^c{a){Ty{xt^) - a)) 

l<ii<...<jfc<n 

f{zi, Zk)pi'^''{npscia){Ty{x[]^') - a), ... , npsc{a){Ty{x[]^) - a))dzi . . . dzk- 

Obviously the same result holds for T!f\x^'^^). And in fact they can hold simultaneously by 
taking the intersection. 

Proof. We will adopt the notation 

;.(T(i)(x(i))) = J2 /(^Psc(«)(T,(x«) - a), . . . ,npM{Ty{xfJ) - a)). 

l<ji<. ..<ifc<n 

Let K^^^ := {z^^^ : {y, z^^^) G Y^'} and i^f ^ := {z^^^ : {y, z^^^) G Y^^}. For any e G (0, 1), 
define the following projections of and Y^^ onto their common domain: 

Yff^ := {y G M^'^'^ : P,[K«] > e} 



F^^; := {y G M^-^" : Pji^f ] > e}, 
^(i)i 



Then for each fixed e, limAr P[Y^^] = 1, which implies there is a sequence e(A^) 1 1 such that 
Y^^'^^ := Yj^l^j^^ has asymptotic probability 1 as well. Therefore, for y G Y^^'^^, 

E,/.nT«(x«)) 

/(^i, . . . ,z,)p(^)(np,,(a)(T,(x«) - a), . . . , np,e(a)(T,(x«) - a)) + o(l) + 0(1 - e(iV)). 

□ 

Thus combined with the result on conditional asymptotic independence, we have for 
a,/3G (-2,2), 

limE,,/„nT(^n^^^^))^?^(T«)2)(a:(2))) 

pi';\ui,...,Uk)f{ ^ A ui -a),...,^— (wfc - a))dui . . . duk 
{[-2,2]fe)t Psc(a) Psc(a) 

pic (^1, • • • , ^0/(^7^(^1 - /3), ■ ■ • , -^T^i^i - (3))dvi . . . dvu 

([-2,2]')^ Psc(p) Psc(/3) 

for y^ G Y^^''^^ This implies unconditional asymptotic independence of the f'^{T^^\x^''^)) 
and (?^(r(2)(x(2))). 

To obtain Theorem 11.11 we need to transfer the independence from indices to the energy 
levels. We will say a pair of energy levels Ei < E2 is covered by a pair of indices Li < L2 
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if the classical locations of [A/,^, Ai^+„+i] and Xi^+n+i] covers Ei and E2 respectively. 
More precisely, let Tl~^ be the inverse cumulative function of the semicircle law on [—2,2], 
then Li covers Ei if 

n-\Li/N) <Ei< n-\{Li +n + 1)/N) 
and similarly for L2 covering E2. Furthermore we say {Ei,E2) is 5-coverd by (Li,L2) if 

Tl Til 

n-\Li/N) + 6-<Ei< n-\{Li + n + 1)/N) - 5-, 

i.e., the covering is padded. 

We will denote by Ts{E) the set of all indices that 5-cover E, and similarly denote 
Zs{Ei, E2) the set of all pairs of indices that 5-cover Ei, E2 respectively. Among Ts{Ei, E2) 
we denote by Qs{Ei,E2) the set of good pairs of indices. Now we have 

Lemma 5.2. Let Ii = [oj — ei,ai + Cj], i = 1,2 (for simplicity), be as in the statement of 
theoremini Denote by K,^s = {(^1,^2) E h x I2 : > e}. Then 

ii,„ ^<^^-'' ^ 



A(/i X I2) 

for every e,6 E (0, 1). Here A denotes Lebesgue measure on M. 

This is proved essentially the same way as the previous lemma. So again we can let e{N) 
be an increasing sequence going to 1, and 5{N) J, 0, with 

lim = 1. 

N A(/i X I2) 

Now to prove Theorem 11.11 it suffices to show that for any 5 > 0, and for asymptotically 
A-almost all {Ei,E2) G /i x I2, the following is true: 
for almost all (Li,L2) G Xs{Ei, E2), 

E[/^,(r(lHa;^'^))^7^,(^f (a;^'^))bLi,yL,+l,l/L,-n,?/L,-n+l] = E/H(Ai, . . . , A^)^H(Ai, . . . , A^) 

where cti = ii22_J^£aI and a2 = "^''"^ ?^-c-2-") _ ^g^^ then let 5 go to 0. But this follows 

^ VLi + l-yLi ' ^ yL2-n+l-yL2-n ^ 

directly from the previous lemma. 
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